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Abstract 


In recent years, the mutual impact of land use change and soil erosion has become a major 
environmental concern. Therefore, this study tried to investigate the changes in different land 
uses and evaluate the impacts of land use changes on soil erosion in Givi city. To achieve the 
goals of the research, first, a land use map was prepared using the object-oriented method for 
the two periods of 2000 and 2021. Then, according to the natural and human conditions of the 
area, other effective factors for erosion in the area were identified and the information layers of 
each of the criteria were prepared by the geographic information system. Evaluation and 
standardization of layers were done using the fuzzy membership function and weighting of the 
criteria was done using the CRITIC method. The final analysis and modeling were done using 
the MABAC method as one of the Multi-Criteria Decision-Making (MCDM) methods. The 
results showed that the largest area in 2000 is related to good and medium pastures, recording 
313.172 and 283.144 square kilometers, respectively. In 2021, the biggest areas are related to 
poor pastures and barren lands, with 335.077 and 329.815 square kilometers, respectively. 
According to the 2000 erosion zoning map, 16.34% and 20.36% of the city and according to the 
2021 erosion zoning, 22.92% and 25.58% of the city’s area are in very high and high erosions, 
respectively. Decreasing good and average pastures and lands with dense vegetation and turning 
them into agricultural areas, poor pastures, and man-made and barren lands have had the 
greatest impact on soil erosion. 
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حغرافیا و مخاطرات محیطی. سال دوازدهم. شمارةٌ جهل و هشتم. زمستان EI‏ صص ۱-2 


مقاله پژوهشی 


بررسی اثرات تغییرات کاربری اراضی بر فرسایش خاک شهرستان گیوی با استفاده از مدل تصمیم‌گیری 
جند معیاره MABAC‏ و تصاویر ماهواره‌ای لندست (سنجنده‌های (TM - OLI‏ 


وب عابدینی - استاد ژئومورفولوژی» دانشگاه محقق اردبیلی» اردبیل» اران 
الناز پیروزی- دکتری ژئومورفولوژی. دانشگاه محقق اردبیلی. اردبیل. oly!‏ 


مهدی فعال نذیری- دانشجوی دکتری ژئومورفولوژی, دانشگاه محقق اردبیلی اردبیل ایران 


تاریخ دریافت: ۱2۰۰/۹/۲۹ تاریخ بازنگری: ۱۶۰۰/۱۱/۳۰ تاریخ تصویب: ۱۶۰۰/۱۲/۱۹ 


حکده 


در سال‌های اخیر, تأثیر متقابل تغییر کاربری اراضی و فرسایش خاک تبدیل به یک نگرانی عمده زیست 
محیطی شده است. ازاین‌رو هدف این مطالعه بررسی تغییرات کاربری‌های مختلف و ارزیابی اثرات 
تغییرات کاربری‌ها بر فرسایش خاک» در شهرستان گیوی است. در راستای دستیابی به اهداف پژوهش: 
ابتدا نقشه کاربری اراضی با استفاده از روش شی گرا برای دو 92 by‏ ۲۰۰۰ و ۲۰۲۱ تهیه شده است. در 
مرحله بعد با توجه به شرایط طبیعی و انسانی محدوده. ple‏ عوامل مؤثر برای فرسایش در منطقه شناسایی 
شد و لایه‌های اطلاعاتی هر یک از معیارهاء توسط سامانه اطلاعات جغرافیایی تهیه گردید. ارزش‌گذاری و 
استانداردسازی لایه‌ها با استفاده از تابع عضویت فازی و وزن‌دهی معیارهاء با استفاده از روش CRITIC‏ 
انجام شد. تحلیل و مدلسازی نهایی با استفاده از روش MABAC‏ به‌عنوان یکی از روش‌های 


تصمیم گیری Ao‏ معیاره (MCDM)‏ انجام ۳ نتایج نشان داد بیشترین ميزان مساحت در Yous his‏ 
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۲ جغرافیا و مخاطرات محیطی شمارة چهارم 


مربوط به مراتع خوب و متوسط به ترتیب با ۳۱۳/۱۷۲ و ۲۸۳/۱۵۶ کیلومترمربع است و در سال ۲۰۲۱ 
بیشترین ميزان Calis‏ مربوط به مراتع ضعیف و اراضی بایر (هموار) به ترتیب با ۳۳۵/۰۷۷ و 4۹/۸10 
کیلومترمربع است. با توجه به نقشه پهنه‌بندی فرسایش سال ۲۰۰۰ به ترتیب ۱۱/۳۶ و ۲۰/۳۹ درصد و با 
توجه پهنه‌بندی فرسایش ۲۰۲۱ به ترتیب ۲۲/۹۲ و ۲۵/۵۸ درصد از مساحت شهرستان در دو طبقۀ 
فرسایش بسیار زیاد و زیاد قرار دارند. کاهش مراتع خوب و متوسط و اراضی با پوشش گیاهی متراکم و 
تبدیل آن‌ها به مناطق کشاورزی» مراتع ضعیف» انسان‌ساخت و اراضی بایر بیشترین میزان تأثیر را بر 
فرسایش خاک داشته است. 


۱- مقدمه 


فرسایش خاک یک مشکل در مقیاس جهانی است که منابع آب‌وخاک را تهدید می‌کند و تغییرات کاربری 
اراضی یکی از عوامل مهم و مؤثر در تشدید خاک است (عابدینی و همکاران. ۱۶۰۲). کاربری اراضی منعکس کننده 
ویژگی‌های تعاملی بین انسان و محیطزیست و توصیف نحوه بهره‌برداری انسان برای یک يا چند هدف بر روی زمین 
است. در حال حاضر تغییر کاربری زمین به‌صورت غیراصولی از مهم‌ترین معضلات اقصی نقاط کشور ایران و جهان 
است و تغییر SI‏ کاربری‌ها WE‏ بی‌برنامه و بدون در نظر گرفتن محدودیت‌های زیست‌محیطی صورت می‌گیرد shes)‏ 
امیر و حاجی قادری. TAT‏ لذاء آگاهی از تغییر و تحولات کاربری اراضی در طول یک دوره زمانی برای برنامه‌ریزان 
Sl pda‏ سار Ge‏ اهمیت است. Bek)‏ و OY SL‏ معمرلا آشکارسازی تخیرات فر اند روط به 
تشخیص مناطق تغییر یافته در دو تصویر و دو زمان متفاوت است (اکسلو و بو" CVV‏ روش‌های متعددی برای 
آشکارسازی تغییرات محیطی به کار گرفته می‌شود که از مهم‌ترین آن‌ها می‌توان به استفاده از تکنیک‌های دورسنجی و 
سیستم اطلاعات جغرافیایی اشاره نمود (جبار و O55‏ ۲۰۰۱؛ وسلس و همکاران؛ ۲۰۰۶: رضایی‌مقدم و همکاران؛ 
۵۹ روش‌های سنتی برای محققین زمانبر و پرهزینه است» ولی تکنیک‌های سنجش‌ازدور با استفاده از تصاویر چند 
زمانه. حجم عظیمی از اطلاعات را در اختیار محقق می‌گذارد (کنس و همکاران“ ۲۰۱۲). همچنین» تصاویر ماهواره‌ای 
به جهت ارائه اطلاعات بهنگام تنوع اشکال» رقومی بودن و امکان پردازش در تهیه نقشه‌های کاربری اراضی از 
اهمیت بالایی برخوردارند. بنابراین می‌توان گفت. استفاده از داده‌های سنجش‌ازدور در محیط GIS‏ می‌تواند شناخت 


مناسبی از چگونگی تغییرات کاربری اراضی ارائه داده و در مدیریت آن راهکارهای مناسبی ارائه دهد (مندوز و 
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سال دوازدهم بررسی اثرات تغییرات کاربری اراضی بر فرسایش خاک .... ۳ 


همکاران» ۲۰۱۱؛ کویین و همکاران» ۲۰۱۶). از سوی دیگر می‌توان فرس‌ایش خاک را یکی از بزرگ‌ترین موانع 
برای دستیابی به توسعه پایدار در کشاورزی و استفاده از منابع طبیعی دانست (سابوتا و همکاران ؛ ۲۰۱۷؛ سینگ و 
FLL‏ ۲۰۱۷). فرسایش خاک فرآیندی است که در آن ذرات خاک سنگدانه‌ها و قطعات سنگی جدا شده و توسط 
آب انتقال می‌یابند و سپس در مناطق جدید قرار می‌گیرند Sul ok)‏ و همکاران؛ ۲۰۱۷). این امر باعث تخریب 
شدید زمین می‌شود و بر باروری خاک بهره‌وری کشاورزی, کیفیت آب و اکوسیستم‌ها و محیط‌های منطقه‌ای تأثیر 
می‌گذارد (آیلو ۲۰۱۵). به‌طور متوسط فرسایش خاک در ایران» سه برابر بیشتر از ple‏ کشورهای آسیایی و ۲۲ برابر 
ميانگین Slee‏ است (مزبانی و همکاران. ۱۶۰۰). بر اساس پژوهش‌های صورت گرفته» کشت بدون عملیات کنترلی 
و همچنین استفادة بدون برنامه‌ریزی از زمین در کنار شسهرسازی» جنگل‌زدایی و کشاورزی غیراصولی از عوامل 
اساسی فرسایش خاک در جهان و به‌ویژه در Ol yl‏ محسوب می‌گردد. پدیدة فرسایش و آثار سوء آن شاید در 
کوتاه‌مدت جندان چشمگیر نباشد. ولی در بلاندمدت بسیار محسوس است. در صورت عدم اقدام» خسارات 
عظیم و جبران‌ناپذیری به همراه خواهد داشت (اوزاحین و ایراغلو ؛ ۲۰۱۹). سیستم اطلاعات جغرافیایی» ابزاری PH‏ 
برای سیستم‌های پشتیبانی تصمیم‌گیری تعاملی برای عملیات مدیریت ریسک هستند (بوخیر و همکاران ۲۰۰) و 
در سال‌های اخیر» جهت مطالعه پدیده فرسایش خاک مورد استفاده فراوان قرار می‌گیرند (آکومار و اراجو ۲۰۱۶ 
آجکومار و همکاران / ۲۰۱۹). به‌علاوه ازآنجایی که مدل‌های مبتنی بر معادله جهانی فرسایش موفق به در نظر گرفتن 
وابستگی متقابل عوامل موثر بر فرسایش خاک نمی‌باشند (نخای" " و همکاران. ۲۰۰۹). لذه برای کشف طیف وسیعی 
از گزینه‌ها از نظر درگیری‌های عینی و معیارهای چندگانه» از روش MCDM‏ استفاده می‌شود (آهر و همکاران ۲ 
۳ در این روش‌های تصمیم‌گیری برای بهینه‌سازی مدل» از معیارهای مختلفی برای افزایش صحت تصمیمات 
استفاده می‌شود (جورجیو و همکاران ‏ ۲۰۱۵؛ عرب‌عامری و همکاران* ۲۰۱۸). لذا می‌توان گفت. با توجه به 


محدودیت‌های روش‌های ستیء که عملتاً وقت‌گیر و هزینه‌بر بوده و اغلب با خطا همراه می‌باشنده استفاده از 
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۴ جغرافیا و مخاطرات محیطی شمارة چهارم 


سیستم‌های اطلاعات جغرافیایی و سیستم‌های تصمیم SS‏ چندمعیاره با یک رویکرد تلفیقی, علاوه بر سودآوری 
می‌تواند باعث تسریع در روند برنامه‌ریزی در تشخیص موارد بحرانی و اضطراری شده و منجر به صدور نتایج 
مناسبی گردد. 
مطالعات فراوانی پیرامون تهیه نقشه‌های کاربری اراضی» پهنه‌بندی فرسایش و نیز نقش کاربری در فرسایش انجام 
شده است. برای مثال؛ اصغری سراسکانرود و همکاران (۱۳۹۷)» اثرات تغیبرات کاربری اراضی را بر روی فرسایش 
g>‏ 85 سیمینه‌رود بوکان با استفاده از روش WLC‏ (طی سال‌های ۲۰۰۰ و ۲۰۱۷ موردبررسی قرار دادند و به این 
نتیجه دست یافتند که بیشترین مساحت این حوضه در iib‏ پرخطر قرار دارد و با افزایش مساحت زراعت eel‏ 
اراضی بای مناطق سنگلاخی و مناطق مسکونی و کاهش مراتع بر میزان فرسایش افزوده می‌شود. میاحی و همکاران 
(Eo)‏ اثر تغییر کاربری اراضی بر فرسایش خاک تالاب هور العظیم را در طی سال‌های ۱۳۹۵ و ۱۳۹۵ موردبررسی 
قرار دادند. نتایج مطالعه OLE‏ داده است در بازه زمانی موردبررسی از مساحت طبقات دارای کلاس فرسایشی کم 
کاسته شده و به مساحت کلاس‌های فرسایشی زیاد افزوده شده است و کاهش پوششگیاهی و مناطق آبی به همراه 
افزایش اراضی بایر نقش مهمی در افزايش درصد فرسایش خاک در محدوده داشته است. عابدینی و همکاران 
(۱۶۰۱) ارتباط شاخحص‌های پوشش گیاهی و ژئومورفیک با مقادیر فرسایش و رسوب در حوزه آبریز کوزه توپراقی 
را موردمطالعه قرار دادند. نتایج ایشان نشان داد که ارتباط میان مقادیر فرسایش و رسوب در زیرحوزه‌های مورد 
مطالعه معنی‌دار است. علاوه براین» مقدار فرسایس و رسوب با مقادیر شیب نیز دارای رابطه معکوس و معنی‌دار 
هستند. هم‌چنین میان شیب و نیز شاخحص قدرت آبراهه رابطه مثبت و معنی‌دار از نظر آماری وجود دارد. سانتوس و 
همکاران" (TY)‏ به بررسی تغییرات کاربری اراضی و نقش آن بر فرسایش در مناطق نیمه‌عشک برزیل» در بازه 
زمانی ۲۰۱۰ تا ۲۰۱۱۱ پرداختند و نتایج به دست آمده نشان داد که سطح اراضی جنگلی و مراتع کاهش aL‏ و ميزان 
فرسایش در این مناطق ۱۰ برابر افزایش داشته است. GAS‏ و همکاران" (۲۰۱۹) تأثیر استفاده از زمین بر فرسایش 
خاک و عملکرد رسوب در اتیوپی. طی سال‌های ۱۹۷۳ تا ۰۲۰۱۵ مورد مطالعه قرار دادند نتایج نشانگر تبدیل سریع 
کاربری‌های با پوششگیاهی طبیعی به کاربری‌های اصلاح شده بشر می‌باشد. کوجو و همکاران" (۲۰۲۰) به بررسی 
تأثیر کاربرق زمین بر فرسایش خاک در کنیا غربی» در سال‌های ۱۹۹۵ و ۲۰۱۷ پرداختند. نتایج نشان داد که ۵۰ درصد 
از تلفات و فرسایش خاک در زمین‌های زراعی» به‌ویژه در مزارعی که در اثر پاک‌سازی و تخریب جنگل‌ها به وجود 
آمده‌انده دیده می‌شود و مناطق دارای پوشش علفی /درختچه‌ای, در رتبه دوم قرار دارند. لی و همکاران" (۲۰۲۱) به 
بررسی تأثیر تغییر کاربری اراضی بر فرسایش رودخانه‌ای در شهرستان Kedong‏ واقع در شمال شرقی چین بین 
Santos et al‏ 1 
Kidane et al‏ 2 


3 Kogo et al 
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سال‌های ۵ تا ۲۰۱۵ پرداختند. نتایج نشان داده است. زمین‌های زیر کشت در بازه زمانی مورد مطالعه به‌طور 
چشم‌گیری افزایش يافته است و کشت زمین در daada‏ به‌ویژه در دامنه‌های بزرگ‌تر از ٤‏ درجه» ممکن است 
سرعت فرسایش رودخانه را تشدید AS‏ 

بخش اعظمی از مساحت شهرستان گیوی» کوهستانی و پر شیب است و نیز اغلب کاربری اراضی بە‌صورت 
زمین‌های کشاورزی» SEL‏ و مراتع می‌باشد که به flo‏ استفاده غیراصولی و بی‌رویه از اراضی» میزان فرسایش 
شهرستان obj‏ می‌باشد و به همین علت از polie‏ زیاد بار معلق رسوبی برخوردار است و این فرآیند در کنار وقوع 
طغیان‌های رودخانه» Col‏ کاهش عمر مفید سد گیوی و کانال‌های واقع بر روی این رودخانه می‌گردد؛ بنابراین» در 
تحقیق حاضر با توجه به اهمیت مطالعه تغییرات کاربری‌ها و نقش آن در فرسایش خاک در گذر زمان و نیز با ASU‏ 
بر بررسی کارایی و دقت نتایج حاصل از به‌کارگیری داده‌های سنجش‌ازدور (RS)‏ سیستم اطلاعات جغرافیایی 
(GIS)‏ و فنون تصمیم‌گیری چندمعیاره (MCDM)‏ با یک رویکرد تلفیقی: تغییرات کاربری اراضی در شهرستان 
گیوی و نقش آن در میزان فرسایش (بین سال‌های ۲۰۰۰ و (VON)‏ مورد مطالعه قرار می‌گیرد. با توجه به این امر که 
تاکنون مطالعه‌ای در زمینه بررسی نقش تغییرات کاربری اراضی بر فرسایش خاک در طی بازه زمانی مورد مطالعه در 
سطح شهرستان گیوی انجام نشده است و به علاوه استفاده از روش ماباک در امر پهنه‌بندی فرسایش خاک برای 
نخستین بار در این مطالعه مورد استفاده قرار گرفته است. نتایج حاصل از این مطالعه می‌تواند برای پژوهشگران, 
مسئولان و سیاست گذاران در راستای پهنه‌بندی فرسایش خاک بسیار جالب باشد و انتظار می‌رود. با شناسایی مناطق 
اولویت‌دار از نظر خطر فرسایش, اقدامات لازم جهت کنترل فرسایش شهرستان گیوی» برنامه‌ریزی بهره‌برداری بهینه 
از خاک و استراتژی‌های مدیریت پایدار زمین به عمل آید. 


۲- مواد و روش‌ها 

۱-۲منطقه مورد مطالعه 

شهرستان گیوی GS)‏ با مرکزیت شهر گیوی به‌طور تقریبی ۱۲۷۹/۷۵ کیلومترمربم» وسعت دارد. این شهر از 
نظر جغرافیایی در tA‏ درجه و ۱ دقیقه تا ÉA‏ درجه و ۲۳ ais‏ طول جغرافیایی و ۳۷ درجه و YV‏ دقیقه تا ۳۷ و OV‏ 
دقیقه عرض جغرافیایی» در فاصله ۸۵ کیلومتری اردبیل قرار دارد (شکل ۱). این شهرستان از شمال با شهرستان 
اردبیل» از شرق با شهرستان خلخال» از غرب و جنوب با شهرستان‌های نير و ale‏ همسایه است و شامل دو بخش 
فیروزآباد و مرکزی و چهار دهستان به نام‌های سنجبد شمالی» سنجبد غربی» سنجبد جنوبی و زرج LT‏ و ۱۲۱ روستا 
ا و ی درا وکر رک کی رای ان که ما E‏ مارت 
دارای اقلیم نیمه حشک خفیف وبا توجه به طبقه‌بندی آمبرژه دارای اقلیم ارتفاعات می‌باشد. ارتفاع متوسط در 


شهرستان گیوی» حدود ۲۳۰۰ متر می‌باشد و بیش‌تر نقاط ارتفاعی آن, در قسمت شمال شرق و جنوب غرب 


۶ جغرافیا و مخاطرات محیطی شمارة چهارم 


قرارگرفته‌اند و از مهم‌ترین ارتفاعات این شهرستان رشته‌کوه تالش يا باغرو داغ می‌باشد. در مناطق کوهستانی 
شهرستان خاک از نوع ریگوسول و لیتوسول بسیار کم عمق تا کم عمق همراه با گراول. ماسه» سیلت رس و قطعات 
زاویه‌دار سنگی متلاشی شده است و در تراس‌ها و حواشی رودخانه‌ها خاک‌ها از لیتوسول نیمه عمیق تا کم عمق 
می‌باشد و در بعضی bw‏ خاک‌ها به‌صورت. نیمه تکامل یافته است. شهرستان گیوی به‌طور عمده از سازندهای 
آتشفشانی و آذرآواری ترشیاری (به‌ویژه ائوسن» الیگوسن و میوسن) تشکیل شده که با رسوب گذاری سنگ‌هایی نظیر 
کنگلومرا» ماسه‌سنگ» مارن و لایه‌های گچ و آهک همراه بوده است. به علاوه» رسوبات کواترنری محدوده؛ شامل؛ 
پادگانه‌های «cd al‏ رسوبات دامنه‌ای و نهشته‌های رودخانه‌ای می‌باشد (شکل ۲ بنابراین» شهرستان کرٹ از مناطق 
فرسایش‌پذیر بوده و انواع فرسایش آبی سطحی» ورقه‌ای. شیاری در آن دیده می‌شود و به همین علت از مقادیر زیاد 


بار Ghee‏ رسوبی برخوردار است. 
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-Y-Y‏ داده‌ها و ابزارهای مورد استفاده 


پژوهش حاضر از نوع کاربردی بوده و روش تحقیق آن تحلیلی مبتنی بر تلفیق آنالیز داده‌هاء سیستم اطلاعات 
جغرافیایی» سنجش ازدور و استفاده از فنون تحلیل چند معیاره است. جهت پردازش تصاویر و تجزیه‌وتحلیل داده‌ها 
نیز از نرم‌افزارهای Excel Idrisi ENVI, Ecognition, Arc GIS‏ بهره‌گیری شده است. جهت بررسی میزان 
فرسایش. در مرحله جمع‌آوری داده‌ها؛ ابتدا باید مؤلفه‌های تأثیرگذار مشخص شده و بر مبنای آن‌ها داده‌ها و اطلاعات 
مورد نیاز گردآوری و طبقه‌بندی شوند. بدین منظون در این مطالعه» ابتدا عوامل موثر (شامل: شیب. لیتولوژی» کاربری 
اراضی؛ بارش فاصله از راه ارتباطی» فاصله از آبراهه و خاک) با توجه به شرایط طبیعی و انسانی منطقه شناسایی 
گردید. در مرحله بعد لایه‌های اطلاعاتی مربوط به هر یک از عوامل. در محیط سیستم اطلاعات جغرافیایی تهیه 
گردید. لایه‌های اطلاعاتی راه‌های ارتباطی و شبکه آبراهه با استفاده از نقشة خحطوط ارتباطی و رودخانه‌های استان 
اردبیل aug‏ گردید. لایه‌های شیب با اسفاده از مدل رقومی ارتفاعی ۱۲/۵ ope‏ دانلود شده از سایت 
(vertex.daac.asf.alaska.edu)‏ تهیه شد. لایه‌های اطلاعاتی مربوط به لیتولوژی (مقاومت سنگ‌ها) ‏ با 
رقومی‌سازی از روی نقشه زمین‌شناسی گیوی» خلخال-رضوانشهر و ماسوله؛ با مقیاس ۱:۱۰۰۰۰۰ تهیه شد. نقشه 
بارش شهرستان go‏ با استفاده از داده‌های ایستگاه‌های هواشناسی و باران‌سنجی و با استفاده از روش درون‌یابی 
۷ ترسیم شد. 


۸ جغرافیا و مخاطرات محیطی شمارة چهارم 


ang -۳-۲‏ نقشه کاربری اراضی 

در تحقیق حاضر با هدف استخراج کاربری اراضی شهرستان GES‏ ابتدا تصاویر ماهواره‌ای لندست از سنجنده‌های 
(TM - OLD‏ برای سال‌های (۲۰۰۰ -۲۰۲۱) در ماه آگوست از سایت زمین‌شناسی آمریکا اخذ شد (جدول۱). برای 
استخراج نقشه کاربری اراضی از ترکیب باندهای مرئی (آبی و سبز و قرمز) و مادون قرمز نزدیک ماهواره لندست. 
جهت طبقه‌بندی استفاده شد. لازم به ذکر است» بر GE‏ طبقه‌بندی پیکسل پایه در طبقه‌بندی شی گرا ترکیب چهار 
باند (۷-1-۳-۲-۱) ماهواره لندست امکان‌پذیر می‌باشد. سپس به‌منظور آماده‌سازی Cp glad‏ تصحیحات هندسی و 
اتمسفری با استفاده از روش ۳1250 روی تصاویر و نرم‌افزار Envi5.3‏ اقدام شد. درنهایت طبقه‌بندی با روش 
co‏ و با الگوریتم نزدیک‌ترین همسایگی توسط نر‌افزار Ecognition‏ صورت گرفت. طبقه‌بندی شیءگرا 
فرایندی است که کلاس‌های پوشش اراضی را به اشیاء تصویری ارتباط می‌دهد این نوع از طبقه‌بندی بر اساس منطق 
فازی استوار است و ارزش عوارض را به ارزش فازی (بین صفر و یک) با درجه عضویت معین برای هر کلاس 
تبدیل می‌کند (رضائی‌مقدم و همکاران, ۱۳۸۹). فرآیند و نتیجه این نوع از تحلیل تصاوین در سه بخش قابل تقسیم 
است که عبارتند ازه سگمنت سازی» طبقه‌بندی فازی و ارزیابی صحت. مبنای طبقه‌بندی شی گرا فرایند قطعه‌بندی 
(سگمنت سازی) می‌باشد. در فرآیند قطعه‌بندی مفسر می‌تواند با در نظر گرفتن پارامترهای شکل, بافت» ضریب 
فشردگی و معیار نرمی» فرایند قطعه‌بندی تصویر را انجام داده و بر اساس نسبت اهمیت هر یک از این عوامل در 
خصوصیات کلاس‌های موردنظر برای استخراج از تصاویر ماهواره‌ای» نسبت اهمیت آن‌ها را در فرایند قطعه‌بندی 
اعمال نماید (فیضی‌زاده و همکاران, ۱۳۸۷). در تحقیق حاضر با تفسیر منحنی‌های انعکاس طیفی و همچنین ارزیابی 
هیستوگرام همبستگی نوارهاء ترکیب نواری ۱-۳-6 به‌عنوان بهترین ترکیب فوق برای طبقه‌بندی انتخاب شد. در 
فرایند سگمنت‌سازی اقدام به اعمال وزن برای ترکیب نواری منتخب گردید و به نوار ۳ (نوار شاخص گیاهی) وزن 
بیشتری از سایر نوارها داده شد و برای این نوار ۶۰ از ۱۰۰ در نظر گرفته شد. در تحقیق حاضر برای قطعه‌بندی از 
روش قطعه‌بندی چند مقیاسه استفاده شد. برای این منظور با تجزیه‌وتحلیل. مقیاس ۸٩‏ و شکل ۰/۱ و ضصریب 
فشردگی ۰/۶ برای قطعه‌بندی انتخاب شد. پارامترهای موثر بر استخراج کلاس‌ها در طبقه‌بندی شی گرا در این 
پژوهش شامل۱- انحراف معیار باندها: انحراف معیار یکی از شاحص‌های پراکندگی است که نشان می‌دهد به‌طور 
میانگین. داده‌ها چه مقدار, از مقدار متوسط فاصله‌دارند در این شاخص از سه پارامتر شامل, باند مادون‌قرمز و قرمز و 
سبز استفاده شد. ۲-شاخحص "NDVI‏ در این پژوهش از شاخحص NDVI‏ به‌منظور جداسازی بهتر پوش شگیاهی از 
بقیه پدیده‌ها استفاده شد.۳-شاخحص PCI‏ ۶سواریانس: واریانس برای محاسبه بافت بر ah‏ ميانگین و پراکندگی در 


1 Normalized Difference Vegetation Index 
2 Precipitation Concentration Index 


اطراف میانگین از ارزش سلول‌ها داخل 'GLCM‏ به کار گرفته می‌شود. ۵-در شاحص هندسی. از پارامترهای مژثری 
چون حداکثر اختلاف. درجه روشنایی. ضخامت» مساحت و فشردگی مورداستفاده قرار گرفت. به‌طورکلی» در 
تعریف شرایط طبقه‌بندی برای هریک از طبقات علاوه بر اطلاعات tb‏ اطلاعات مربوط به موقعیت. محتوا؛ بافت 
و ple‏ ویژگی‌های بصری نیز استفاده شده است به‌عنوان نمونه بافت منحصربه‌فرد محدوده مورد مطالعه برروی 
تصویر شناسایی bol‏ را از سایر طبقات کاربری اراضی میسر می‌ساخت. بنابراین در طبقه‌بندی شیءگرای تصویر 
برای این عامل ارزش بیشتری داده شد. همچنین مشخص شد که اراضی تحت آیش کشاورزی دارای اطلاعات طیفی 
یکسانی با طبقه کاربری» اراضی pb‏ هستند و با طبقه‌بندی بر اساس ارزش عددی پیکسل‌هاء تفکیک این دو طبقه 
کاربری اراضی از هم به‌حوبی امکان‌پذیر نیست» اما با توجه به اینکه اراضی تحت آیش کشاورزی دارای اشکال تقریاً 
منظم هندسی هستند. تنها شکل چندضاعی آنها تفکیک این دو طبقه را ممکن می‌سازد. لذا عامل شکل, به‌عنوان 
عاملی تأثیرگذار در تفکیک طبقه اراضی تحت آیش از اراضی بای شناخته شد که در تعریف شرایط طبقه‌بندی در 
تفکیک این دو طبقه از هم استفاده شد. برای سایر طبقات کاربری اراضی نیز متناسب با ویژگی‌های بصریشان علاوه 
بر اطلاعات طیفی از ple‏ پارامترها نظیر شکل» بافت. تن رنگ و زمینه استفاده شد. بررسی نتایج طبقه‌بندی یکی از 
مراحل مهم ULC ol‏ ارلئه نتایج طبقه‌بندی بدون هیچ‌گونه پارامتری که کیفیت با صحت این نتایج را بیان کند. از 
ارزش آن‌ها می‌کاهد و در بعضی مواقع آن‌ها را بدون استفاده می‌کند. دو پارامتر اصلی که برای ارزیابی صحت 
طبقه‌بندی مورد استفاده قرار گرفته است. صحت کلی و ضریب LIS‏ می‌باشند. 


حدول \- خصوصیات تصو پر سنحنده (TM-OLI)‏ از ماهواره لند ست (منبع: متادیتای تصاویر) 


zj‏ تاريخ اخذ ردیف /گذر | تعداد باند قدرت تفکیک زاویه ازیموت خورشید زاویه ارتفاع خورشید 
لندست YN ۱ TM‏ 1 باند ۱۱۳/۰ 04/46 
۱/۳۶ ۰ متر 
لندست Yee ۱/۶ OLI‏ ۷ باند \V/EY‏ 0۹/4 


Jol» -٤-۲‏ روش Any lie) MABAC‏ جند شاخصة محدوده تقریبی مرزها) 
۱- تعیین معیارها و گزینه‌های پژوهش. ۲- ماتریس تصمیم در این روش به‌صورت معیار-گزینه است؛ یعنی یک 
ماتریس که ستون‌های آن را معیارهای مسئله و سطرها را گزینه‌ها تشکیل می‌دهند و هر سلول نیز در واقع امتیاز هر 


mV)‏ تعداد کل گزینه‌های جایگزین را نشان می‌دهد. پس از محاسبه مقدار gi‏ ماتریس محدوده‌های تقریبی مرزی 


1 Gray Level Co Gray Level Co-Occurrence Matrix 
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با توجه به معیارهای G‏ با فرمت ۷*1[ شکل می‌گیرد و 1 تعداد کل معیارهایی را که گزینه‌های ارائه شده برای آنها 
انتخاب می‌شود را نشان می‌دهد؛ cyl ply‏ اگر ٩‏ معیار داشته باشید یک ماتریس 61×1 به‌صورت زیر خواهیم داشت 
(رابطه ۲): 

9: = (IT, w” ۱ رابطه‎ 

G = [g1.g2.....gn] ۲ رابطه‎ 

-محاسبه فاصله گزینه‌ها تا موز ناحیه شباهت: در این بخش با استفاده از رابطه ۳ فاصله گزینه‌ها را تا ناحیه g‏ 
به دست می‌آید. در واقع 4b‏ ماتریس وزندار را از ماتریس 2 کم کرد. 

Q=V-G all, 

پس از مشخص شدن ماتریس Q‏ می‌توان با استفاده از حد بالایی مساحت (GH)‏ و حد پایین مساحت (-6) 
وضعیت هر گزینه را مشخص کرد. بر این اساس گزینه AI‏ متعلق به اجتماع مجموعه مذکور اسست. حد بالای 
مساحت (G+)‏ ناحیه‌ای است که گزینه ou!‏ آل مثبت در آن قرار دارد و حد پایین مساحت (G-)‏ ناحیه‌ای است که 
گزینه ضد ایده‌آل در آن قرار دارد. میزان تعلق گزینه Al‏ به اجتماع بالا بر اساس رابطه 6 به دست می‌آید. برای 
انتخاب گزینه ۸به‌عنوان بهترین فرم از مجموعه. لازم اسست که حداکثر معیارهای ممکن به منطقه تقریبی فوقانی 
(G+)‏ تعلق داشته باشد. مقدار بالاتر +166 نشان می‌دهد که گزینه جایگزین به گزینه جایگزین ایده آل نزدیک‌تر 
است. درحالی‌که مقدار پایین‌تر GiEG—‏ نشان می‌دهد که گزینه جایگزین» به گزینه cp Sle‏ ضد ایده آل نزدیک‌تر 


است. 


Gt if qij >o 
۸ 6106 if Wj=o ۶ abl, 
GT if Qij <o 


۷- انتخاب گزینه بهینه: در روش ماباک (MABAC)‏ با استفاده از رابطه ۵ امتیاز نهایی هر گزینه را مشخص 
کرده و بر اساس آن گزینه‌ها رتبه‌بندی می‌شوند. 
رابطه ۵ Si = J(qij);ji = 1.2 ....njj = 1.2....m‏ 
محاسبه polis‏ توابع معیار» توسط گزینه‌ها به‌عنوان مجموع فاصله‌های جایگزین از مناطق تقریبی مرزی Gi‏ به 
دست می‌آید. با جمع کردن pole‏ ماتریس 0) در هر سطر مقادیر نهایی تابع معیار گزینه‌ها به دست می‌آید. در این 
رابطه n‏ تعداد معیارها را نشان می‌دهد و 10 تعداد گزینه است (علی‌نزاد و خلیلی ۰ ۲۰۱۹؛ وانگ و همکاران ؛ ۲۰۲۰). 


1 Alinezhad & Khalili 
2 Wang et al 


0-7- وزن‌دهی CRITIC‏ (اهمیت معیارهاء از طریق همبستگی بین معیارها) 

با توجه به اینکه در وزن‌دهی CRITIC‏ وزن معیارهاء با توجه به نقشه‌های معیار و میزان همبستگی» تضاد و 
انحراف gles‏ ين لایه‌های اطلاعاتی مورد استفاده تعیین می‌گردد استفاده از این روش در وزن‌دهی معیارها در 
پژوهش حاضر می‌تواند معضل استقلال صفات از یکدیگر که به هنگام مقايسة زوجی در چارچوب روش فرایند 
E ies‏ یل کان byw‏ کا کر اب وان بیان ا اسراف 
معیار معیارهای موردبررسی» ماتریس متقارنی به ابعاد mxm‏ ایجاد می‌گردد که شامل ضرایب همبستگی بین بردارهای 
فشک as‏ انست, با تین پازآترهای E‏ تین مار S VF‏ فیکر از رو رانظا ابد 
می‌شود. 

Cjk = Xk=1(1 - rjk) ٦ رابطه‎ 

که در آن Cik‏ معرف مجموع تضاد sles‏ با معیارهای >1 است که از 1 = ) شروع شده و تا Mm‏ <1 ادامه دارد و 
زه همبستگی بین دو معیار و زرا نشان می‌دهد. میزان اطلاعات عامل زرا می‌توان با استفاده از رابطة V‏ محاسبه 
نمود. 

Cj <<, )1 - rjk) ۷ رابطه‎ 

که در Cj ul‏ معرف میزان اطّلاعات معیار ز و انحراف معیار در مقادیر مربوط به عامل یا معیار زرا OLE‏ می‌دهد. 
با توجه به روابط یادشده معیارهایی که دارای Cj‏ بیشتری باشند وزن زیادی به خود احتصاص خواهند داد. وزن هر 
fale‏ مانند از رابطةٌ A‏ تعیین می‌گردد. 


Cj 
= ۸ ابطه‎ 
از‎ a Ck a 


که در آن Bre Wj‏ وزن معیار زو SLE CK‏ میزان اطلاعات مجموع معیارهای k‏ است که از 1 = k‏ شروع 


شده و تا 110 ادامه دارد. لذا؛ وزن نهایی هر معیارء از تقسیم میزان اطلاعات هر معیار بر مجموع میزان کل اطلاعات 
تمامی معیارها به دست می‌آید (علی‌نژاد خلیلی. ۲۰۱۹؛ تاش و آیتک آدالی؛ ۲۰۱۹). در راستای انجام مراحل وزن 
دهی کرتیک ابندا در محیط ادریسی انحراف معیار و همبستگی Ole‏ لایه‌های اطلاعاتی مربوطه به هر یک از عوامل با 
استفاده از تابع REGRESS‏ به دست آمد و در مرحله بعد» با توجه به روابط مربوطه نسبت به محاسبه ماتریس 
تضاد. میزان اطلاعات و وزن نهایی در اکسل اقدام گردید. 

در پیاده‌سازی تمامی الگوریتم‌های مورد استفاده در پژوهش حاضر و جهت دستیابی به الگوی پهنه‌بندی مکانی 
فرسایش در حوضه گیوی‌چای, مرحله استانداردسازی, با توجه به تابع عضویت فازی انجام شده است (جدول AY‏ 


استاندارد نمودن داده‌ها به معنی همسان‌کردن دامنة تغییرات داده‌ها بین صفر و یک و یا یک دامن مشخص دیگر است 


1 Tus & Aytaç Adalı 


1 جغرافیا و مخاطرات محیطی شمارة چهارم 


در مجموعه‌های فازی» بیشترین ارزش یعنی مقدار یک به حداکثر عضویت و کمترین ارزش یعنی صفر به حداقل 
عضویت در مجموعه تعلق می‌گیرد (بایوم‌گرتل و همکاران؛ ۹ حسینی و همکاران " ۰۹ فلوچارت مربوط 


به فر آیند پژوهش, در شکل ۲ نمایش داده شده است. 


جدول ۲- نوع als‏ فازی‌ساز لایه‌های موضوعی تأثیرگذار بر فرسایش خاک شهرستان گیوی 


i 


تهیه نقشه خروجی با روش MABAC‏ 
وزن دهی با روش CRITIC‏ 


ارزش گذاری و استانداردسازی با استفاده از 


E 


پوش گیاهی متراکم زمین هموار | مناطق آبی | 


شکل ۳- فلوچارت مراحل تحقية 


نتشه کاربری اراضی (۲۰۳۱-۲۰۰۰) | مرتع 


زراعت آبی 


متغیر توضیحات نوع تابع فازی 
شیب کلاسه‌بندی مجدد و اختصاص کدها ر اساس اهمیت آن در وقوع خطر فرسایش و سپس فازی a‏ بزرگ فاری 
کردن. 

کاربری اراضی اختصاص کدها بر اساس اهمیت آن در وقوع خطر فرسایش و سپس فازی کردن. تابع بزرگ فازی 
فاصله از راه ارتباطی در مناطق مجاور راه ارتباطی» پتانسیل وقوع مخاطره فرسایش افزایش می‌یابد. تابع نزدیک فازی 
لیتولوژی اختصاص کدها بر اساس اهمیت آن در وقوع خطر ز فرسایش و سپس فازی کردن. تابع بزرگ فازی 
خاک اختصاص کدها بر اساس اهمیت آن در وقوع خطر فرسایش و سپس فازی کردن. تابع بزرگ فازی 

فاصله از رودخانه در مناطق مجاور رودخانه پتانسیل وقوع مخاطره فرسایش افزایش می‌یابد. تابع نزدیک فازی 
بارش با افزایش مقدار بارش. خطر فرسایش افزایش می‌یابد. تابع حطی افزایشی 


زراعت دیم 


1 Baumgertel et al 
2 Hoseini et al 


سال دوازدهم بررسی اثرات تغییرات کاربری اراضی بر فرسایش خاک .... ۳ 


۳- نتایج و بحث 


دا قهه تفه pal) ea‏ 

کاربری‌های موردنظر با پوشش‌اراضی آن‌ها در ده کلاس؛ شامل: مناطق مسکونی آب (رودخانه‌ها و دریاچه‌ها» 
زراعت ديم زراعت cel‏ پوش شگیاهی متراکم (اراضی با پوشش درخت طبیعی. مانند: جنگل و مصنوعی مانند: 
باغات احداث شده گردو» به و ... )» اراضی کوهستانی (پهنه سنگی و صخره)» مراتع حوب. مراتع متوسط مراتع 
ضعیف و اراضی بایر هموار (اراضی SE Col‏ لخت بدون پوشش گیاهی). طبقه‌بندی شدند. در این پژوهش به‌منظور 
بررسی صحت طبقه‌بندی, از نرم‌افزار 5.3 Envi‏ استفاده شد. لذا جهت محقق نمودن این عمل از روش پیکسل پایه 
و الگوریتم بیشترین فراوانی بهره گرفته شد. طبقه‌بندی الگوریتم بیشترین فراوانی دقیق‌ترین و پراستفاده‌ترین روش 
نظارت‌شده است. فرآیند طبقه‌بندی فرض را بر این می‌گذارد که داده‌های آماری آموزشی برای هر کلاس در هر باند. 
به‌صورت نرمال توزیع شده‌اند؛ بنابراین برای داده‌های آموزشی با هیستوگرام‌های n‏ مدی در یک UL‏ منفرد مناسب 
نیستند در چنین شرایطی هر یک از مدها احتمالاً نماینده کلاس‌های جداگانه‌ای هستند و لازم اسست نرم‌افزار 
طبقه‌بندی کننده» برای آن‌ها مجدداً آموزش داده شوند. بهاین‌ترتیب» هیستوگرام‌های تک مدی با توزیع گوسن که 
شرط توزیع نرمال را برآورد می‌کند به وجود می‌لید. بدین صورت که ابتدا برای هر کاربری به تعداد 1۵ نمونه 
(به‌طور میانگین) از سطح محدوده برداشت شد و سپس در مرحله بعد با طبقه‌بندی شی گرا و بیشترین فراوانی» تک 
به تک کاربری‌ها با یکدیگر مقایسه و بررسی گردید. نتایج جدول ۳ نشان می‌دهد که در سال ۲۰۰۰ صحت کلی ۸٩‏ 
درصد و ضریب LE‏ به‌طور میانگین برابر با ۰/٩۱‏ می‌باشد و در سال ۲۰۲۱ صحت کلی به‌طور میانگین ٩۳‏ درصد و 
ضریب US‏ به‌طور میانگین برابر ۰/۸۷ می‌باشد و با توجه به صحت IS‏ و ضریب LIS‏ نقشه‌ها و با استناد و مقایسه با 
منابعی نظیر (لایل‌سند ۲۰۰۸؛ کنگالتن » ۲۰۰۹) که ضرایب صحت IS‏ بزرگ‌تر از ۸۵ درصد را قابل‌قبول عنوان 
etles S‏ نتایج حاصله از طبقه‌بندی کاربری‌ها چه از نظر تک‌تک کاربری‌ها و چه از نظر مجموع صحت و آمار کاپا از 
صحت bb‏ قبولی در ارتباط با اطلاعات تولید شده برخوردار شده است. درنهایت لایه‌های به دست آمده برای 


محاسبه مساحت کاربری‌ها و تهیه نقشه خروجی مناسب به نرم‌افزار «Arc GIS‏ انتقال داده شدند. 


1 Lillesand 
2 Congalton 


۱۴ جغرافیا و مخاطرات محبطی شمارة چهارم 


جدول ۳- بررسی صحت طبقه‌بندی کاربری اراضی 


ضریب us‏ صحت کلی ماه سال 
Aa% ۰۹۱‏ آگوست Yous‏ 
ay ۰/۸۷‏ آکوست yey)‏ 


E E atlas pale al E‏ تالف کر وال 10۸۷ مرو باه ری 

pial) ر نقشسه کاربری‎ ERB) متوسط به فرئیب با ۳۱۳/۱۷۹ و ۲۸۳/۱۵6 کیلومترمريم می‌باشاد‎ rg ly 
۳۷۹ /۸۱۵ به ترتیب کاربری‌های مراتم ضسعیف و زمین‌های پست و هموار با مساحت ۳۳۵/۰۷۷ و‎ ۱ 
کیلومترمربع» بیشترین مساحت را به خود اختصاص داده‌اند (شکل ۵). کاربری‌های مناطق آبی و انسان‌ساخت کمترین‎ 
مقدار مساحت را در هر دو دوره شامل می‌شوند. با توجه به نتایج مطالعه. شاهد کاهش مراتع متوسط مراتع خوب»‎ 
شاهد افزایش کاربری‌های مراتع ضعیفه زراعت‎ Ss پوشش گیاهی متراکم و زراعت آبی از یک طرف و از طرف‎ 


۴۸۱۰۱۰۲۲ 


نقشه کاربری اراضی شهرستان گیوی 
سال Yess‏ 


> مرتع متوسط‎ 
D انان ساخت‎ ge ا‎ sway 
D متراکہ‎ oe He : 


شکل -٤‏ نقشه کاربری اراضی سال ۲۰۰۰ شهرستان گیوی 


سال دوازدهم 


۴۸۰۴۰۰ ۴۸۰۵۰۰ 


نقشه کاربری اراضی شهرستان گیوی 


< مناطق آبی‎ S متراکم‎ ss ve 
D مرتع ضعیف 3 مناطق کوستانی‎ 


D at زراعت‎ 


۴۸۰۲۰۱۰۲۲ ۴۸۰۳۰۰۳۲ 


FAN "E 


شکل ۵- نقشه کاربری اراضی سال ۲۰۲۱ شهرستان گیوی 


FA! "E 
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جدول ۶- مساحت و میزان تغییرات هر نوع از کاربرهای اراضی در شهرستان گیوی در سال‌های ۲۰۰۰ تا ۲۰۲۱ 


نوع کاربری 
مراتع ضعیف 
مراتع متوسط 
مراتع خوب 
زراعت آبی 
زراعت دیم 
زمین هموار 
اراضی کوهستانی 
پوشش گیاهی متراکم 
مناطق انسان‌ساحت 


مناطق آبی 


در سال ۲۰۰۰ مساحت به 
کیلومترمربع 
1/۹4۳ 
YAY/\ ££‏ 
۱۳۱/۷۲ 
LATAK‏ 
۱۱ 
00/4 
۱۹/۳۹۰ 
۱۳۳/۹۳۰ 
WEAN‏ 


۱۳/۸6۸ 


در سال ۲۰۲۱ مساحت به 
کیلومترمربع 
۳۳۵/۷ 
VENA‏ 
YAVEVO‏ 
۱۳/۳۹۰ 
۱-۳ 
3/۸10 
very‏ 
40/£AY‏ 
V4‏ 


۳/۳۷۵ 


میزان تغییرات به 
کیلومترمریع 
در 
\TE/TAE‏ 
yay‏ 
YN 000‏ 
ADI‏ 
VET \‏ 
۱۳ 
\TANEEV‏ 
۱۳/۳۹ 


ONY 


روند تغییرات 


افزایشی 


۶ جغرافیا و مخاطرات محیطی شمارة چهارم 
۲-۳- پهنه‌بندی فرسایش خاک 
جهت پهنه‌بندی فرسایش با توجه به نقشه کاربری مربوط به دو دوره زمانی و نیز نقشه ساير معیارهای مطرح 
(شکل CV‏ ابتدا وزن‌دهی معیارها انجام گردید. مفروضات پایه‌ای (ماتریس همبستگی ماتریس تضاد. انحراف معیارء 
میزان اطلاعات) و همچنین» وزن نهایی حاصل از وزن‌دهی کرتیک در بین معیارهای مطرح در پهنه‌بندی فرس‌ایش 
خاک مورد مطالعه در جدول ۵. آورده شده است و با توجه به نتایج حاصله به ترتیب معیارهای شیب کاربری 
ارضی و خاک در هر دو دوره مورد مطالعه. بیشترین وزن و اهمیت را دریافت کردند. سپس با انجام مراحل روش 
MABAC‏ نقشه‌ی پهنه‌بندی فرسایش شهرستان گیوی, در ۵ طبقه با فرسایش بسیار زیاد تا فرسایش بسیار کم در 
دوبازه زمانی ۲۰۰۰ و ۲۰۲۱ به دست آمد (شکل ۷ و شکل (A‏ مساحت و تعداد پیکسل‌های مربوط به هر یک از 
طبقات خطر نیز در جدول A‏ ذکر شده است. 


۴۸۲ FAY" "E faero "E FACA! "E 


z 

z G 
$ 2 + Z 
å] > care و‎ 
> é = % ai 2 
£ راهنمای نقشه‎ + = 

محدوده مورد مطالعه > 

EEST 

Do بسیار نامقاوم‎ 
3 Drala 3 
9 3 R z ] |۷۵ E 
3 SD مقاومت متوسط‎ žė ? 
È ge Meer | 
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مقاوم D‏ 
بسيار مقاوم BD‏ , 


E o-i 


کیلومتر | 
۳۰ 


FAFA "E‏ ۴۸۹۳۵ دار را دای 


FA HE FATT ME FAT "E 


۱۳۷۰۵۰۰۸۸ 


FY." 


very! "N 


FA "E FAY "E FACTO! "E FA E FACT "E ۴۸۹۳۵۱۰۷ ع‎ 
نقشه بارش نقشه خاک‎ 
= شهرستان گیوی‎ 
2 
2 - n? = 
S tas K 
i - ۰ Pas — es 
را ای نقشه‎ ane ای‎ r را‎ 
FD محدوده مورد مطالعه‎ FD محدوده مورد مطالعه‎ 
میزان بارندگی ارزش‎ 
z E os 
> as برحسب میلی متر‎ ۳ 
z] شترین :۶۱۸ چ اینسپتی سول‎ 
z D رخنمون سنگی-‎ re = دیشر‎ z 
کمترین:۳۲» تا = انتی سول‎ | 
اینسپتی سول‎ ۲ 
ad — 
1 کیلومتر | —— کیلومتر‎ 
۷ ۲۵.۵ FF ۱۴ ۲ ۸ 
دی‎ pa ry "E fa ra "E F۸44 "E Farry "E ۴۸۳۳۵۰ 


شکل *- نقشه معیارهای مطرح در پهنه‌بندی فرسایش شهرستان SIE‏ 


جدول ۵- مجموع تضاد. انحراف معیار. میزان اطلاعات و وزن نهایی معیارهای مطرح در پهنه‌بندی فرسایش 


شهرستان گیوی 


2 معیار مجموع تضاد انحراف معیار میزان اطلاعات وزن نهایی 
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راهنمای نقشه 


قرسایش توس ie‏ موده CD‏ 


فرسایش زیاد اروش 


is فرسایش خاک سال‎ E سار کہ ارم‎ ed E 
aP فرسایش كم‎ 
شهرستان گیوی‎ 


FARO OPE 


راهنمای نقشه 
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فرسایش زياد ارزش 
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۲۰۲۱ فرسایش خاک سال‎ or ° y فرسایش کم‎ 
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شکل ۸- نقشه پهنه‌بندی فرسایش سال ۲۰۲۱ شهرستان گیوی با استفاده از روش MABAC‏ 


حدول -I‏ اطلاعات طبقات فرسایش خاک در سال‌های ۲۰۰۰ و ۲۰۲۱ شهرستان گیوی 


طبقة فرسایش بسیار زیاد زیاد متوسط کم بسیار کم 

Wy CLO YAV/\YY ۱۳/۵۸۷ 1۹7160 سال ۲۰۰۰ مساحت به کیلومترمریع‎ 
Vik YE/AY YEE ۳/۳۹ \ VTE مساحت به درصد‎ 

MLA ۱۳۱۳۸/۳۳۳ YEVVYA ۱۳/۳۹ ۱۳/۳۶ مساحت به کیلومترمربع‎ Sadie 
voy \V/AE 17/۹ YO/OA YY/4Y مساحت به درصد‎ 


می‌توان گفت که در نقشه فرسایش هر دو دورهٌ مطالعاتی. به‌طور عمده مناطق با فرسایش بسیار زیاد و زیاد در 
کاربری‌های زراعی» مراتع ضصعیف مناطق انسان‌ساخت و اراضی 2b‏ قرار دارند. با توجه به نقشه‌های پهنه‌بندی 
فرسایش در محدوده مطالعاتی» در سال ۲۰۰۰ مساحت طبقه فرسایش بسیار زیاد و زیاد. ۲۰۹/۰3۵ و ۲2۰/۵۸۷ 
کیلومترمریع بوده اسست که مقدار این طبقات خطر در سال ۲۰۲۱ به ترتیب؛ ۲۹۳/۳۵۶ و ۳۲۷/۳6۲ کیلومترمربع 
افزایش یافته است. علت افزایش مناطق دارای احتمال فرسایش را می‌توان در کاهش مراتع حوب مراتع متوسط و 
اراضی با پوششگیاهی متراکم و تبدیل آن‌ها به مناطق کشاورزی, انسان‌ساخت. مراتع ضعیف و اراضی بایر دانست. 
اصغری سراسکانرود و همکاران (۱۳۹۷» اصغری سراسکانرود و همکاران (۱۳۹۸» میاحی و همکاران (۱۶۰۰)» 
سانتوس و همکاران OLS (THY)‏ و همکاران (۲۰۱۹)» کوجو و همکاران (۲۰۳۰)» لی و همکاران (۲۰۲۱» نیز 
افزایش مساحت اراضی کشاورزی. اراضی pl‏ و مناطق مسکونی و در مقابل کاهش مراتع و پوششگیاهی را از 
دلایل اصلی افزایش پتانسیل فرسایش خاک مطرح کرد‌اند. 

قابل‌ذکر می‌باشد که بیشتر اراضی مرتعی توسط روستایبان منطقه به اراضی کشاورزی و به‌عصوص به 
دیم‌زارهای کم‌بازده. تغییر کاربری داده‌اند. چون اکن اراضی کشاورزی منطقه بر روی دامنه‌های شیب‌دار قرار 
گرفته‌اند و بر اثر تغییر کاربری از اراضی جنگلی و مرتعی به کشاورزی تبدیل يافته‌اند. به مرور زمان, خاک و 
مواد حاصل خیز این اراضی بر اثر فرسایش از دسترس خارج می‌شود. WW‏ با توجه به کوهستانی بودن Ailes‏ مورد 
مطالعه. با تبدیل مراتع اطراف روستاها به دیم‌زارها و شخم خوردن مراتع در شیب‌های تند. مواد مغذی درون خاک 
شسته شده و قدرت حاصلخیزی خاک پس از چند سال از دست رفته و احتمال فرسایش بیشتر می‌شود. همچنین؛ 
تغییر کاربری اراضی مرتعی به‌ویژه در مناطق پر شیب کوهستانی علاوه بر فرسایش خاک به کاهش میزان تغذیه 
آب‌های زیرزمینی؛ بروز سیلاب‌های مخرب نیز منجر خواهد شد و کیفیت پویای خاک را تحت تأثیر قرار می‌دهد. 
بايد توجه داشت که پوششگیاهی طبیعی» با نگهداشت قطرات باران باعث کاهش انرژی جنبشی باران شده که باعث 
افزایش پایداری خاک شده و در مقابل فرسایش نیز محافظت می‌شود که با تغییر کاربری موردنظر چنین مکانیسمی از 
بین می‌رود. همچنین» علاوه بحث تغییر کاربری مراتع به زمین‌های زراعی» موضوع دیگری که در رابطه با کاهش 


مساحت ele‏ خوب و متوسط در شهرستان گیوی مشهود می‌باشد؛ بحث تبدیل ely‏ خوب و متوسط به مراتع 
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ضعیف می‌باشد. به علت اينکه واردات خوراک دام هرساله با افزایش قیمت. مواجه می‌باشد. دامداران برای آنکه در 
هزینه‌های خوراک دام صرفه‌جویی نماینده اقدام به چرای بی‌رویه دام‌ها در مراتع می‌نمایند لذا بخشی از خسارت 
وارد شده به SE‏ و پوششگیاهی مراتع شهرستان گیوی, ناشی از عدم رعایت زمان مناسب چرای دام ede‏ توجه به 
ظرفیت چرا و غفلت از اثرات چرای سنگین بر این اراضی بوده است. 

البته بخشی از نابودی پوششگیاهی محدوده در اثر احداث سد گیوی و تخریب پوشش گیاهی طبیعی منطقه. اعم 
از باغات و مراتع بوده است. سد مخزنی گیوی به‌منظور تنظیم آب رودخانه با هدف کنترل سیلاب و نیز تأمین آب 
شرب صنعت و بهبود اراضی منطقه» با ۵۷ میلیون مترمکعب حجم کل مخزنی و ۷۱ میلیون مترمکعب. حجم تنظیمی 
بر روی رودخانة گیوی‌چای. احداث شده است. این سد از نوع سنگریزه‌ای با iia‏ رسی است و دو تونل انحراف و 
سرریز از نوع آزاد دارد. ه‌علاوه سطح دریاچۀ سد ۶۰ هکتار تراز تاج سد ۱۵۵۲ متر و تراز نرمال OT‏ ۱۵۶۸ متر است 
(سازمان OT‏ منطقه‌ای اردبیل» ۱۳۹۲). کاهش میزان مراتع» پوشش گیاهی متراکم و اراضی زراعی در بالا دست سد و 
درنتیجه از بین رفتن پوشش طبیعی منطقه. تغییرات بسیار زیادی را در سیستم‌های مورفولوژیک منطقه ایجاد می‌نماید 
و این تغییرات عمدتاً در جهت افزایش میزان فرسایش و رسوب می‌باشد. به علاوه باید در نظر داشت که افزایش 
مساحت اراضی بایر که به دنبال تخریب اراضی با پوششگیاهی متراکم» مراتع و زراعت آبی در شهرستان رخ داده 
است. به همراه تبدیل اراضی با پوشش گیاهی متراکم و مراتع. به مناطق مسکونی نیز به افزایش سطوح با فرسایش 
زیاد و بسیار زیاد. از نظر فرسایش خاک در سطح شهرستان منجر شده است. 

با توجه به نقشه پهنه‌بندی فرسایش خاک شهرستان گیوی. می‌توان گفت؛ وضعیت توپوگرافی و شیب یکی از 
مهم‌ترین عوامل دخیل در فرسایش خاک محدوده محسوب می‌گردد. معمولاً بر روی دامنه‌های با شیب بیش از ۶۰ 
درصد ضخامت خاک در محدوده کاهش می‌بابد و نقش کمتری در تولید رسوب ناشی از وقوع فرسایش دارند و 
مناطق با فرسایش زیاد و بسیار زیاد. به‌طور عمده در شیب ۱۵ تا ۶۰ درصد قرار دارند. به علاوه در مناطق کوهستانی 
شهرستان خاک از نوع ریگوسول و لیتوسول بسیار کم عمق تا کم عمق همراه با گراول. ماسه. سیلت. رس و قطعات 
زاویه‌دار سنگی متلاشی شده است و در تراس‌ها و حواشی رودخانه‌ها خاک‌ها از لیتوسول نیمه عمیق تا کم عمق و 
مستعد فرسایش می‌باشد و در بعضی نقاط خاک‌ها به‌صورت. نیمه تکامل aL‏ است. از دیگر معیارهای مهم که در 
ارزیابی فرسایش خاک موردبررسی قرار می‌گیرد لیتولوژی می‌باشد. در شهرستان گیوی, واحدهای مختلف 
سنگ‌شناسی منطقه به‌صورت آتشفشانی و رسوبی گسترده شده است. ازجمله این واحدها می‌توان به گدازه‌های 
آندزیتی, گدازه‌ها و گنبدهای داسیتی» تراکی آندزیتی, بازالت و مواد آتشفشانی ایگنمبریت و توف به همراه سنگ‌های 
رسوبی تراورتن. پادگانه‌های آبرفتی جوان و قدیمی که از تخریب و فرسایش سنگ‌های آتشفشانی حاصل شده‌اند. 


اشاره کرد. هر یک از این سازندها از مقاومت متفاوتی برای فرسایش برخوردارند. بررسی نقشه نهایی و مقايسهٌ آن با 


عامل لیتولوژی OLE‏ می‌دهد. سازندهای نامقاوم, که عمدتاً سازندهای سست کواترنری هستند (Qal Qt: Qtr)‏ و 
مارن‌های خاکستری» رس و قرمز ژیپس‌دار (Nems)‏ که بخش اعظم منطقه را پوشانده‌اند. دارای حساسیت بالایی 
نسبت به فرسایش می‌باشند. مشخصات بارندگی مانند شدت. مدت زمان و حجم در ایجاد فرسایش خاک نقش 
مهم و اساسی دارد. مقدار بارش سالیانه در محدوده مطالعاتی ۲۶۳-۷۲۰ میلی‌متر و میانگین بارندگی سالانه حدود 
۰ میلی‌متر بوده و بیشترین مقدار بارش شهرستان مربوط به ماه‌های اردیبهشت و فروردین می‌باشد و Gols‏ با 
بارش زیاد عمدتاً در نواحی کوهستانی و مرتفع گسترده شده‌اند. مناطق با پتانسیل بالای فرسایش, به‌طور عمده در 
مناطق با بارش (۳۰۰-۵۰ میلی‌متر) قرار دارند. همچنین» سطح زیادی از اراضی با فرسایش زیاد و بسیار زیاد در 
حواشی رودنحانه‌های اصلی موجود در شسهرستان (گیوی‌چای, گزاز چای, آرپاچای» نیلق چای) قرار دادند. 
رودخانه‌های شویرچای» سنگورچای. گیوی‌چای نیلق‌چای, لیکوان‌چای از رودخانه‌های دائمی و پرآب شهرستان 
گیوی می‌باشند که شعبات رودخانه‌ها به‌راحتی می‌تواننده مواد سطحی و خاک را شسته و با خود حمل نمایند و با 
فرسایش زیاد و بسیار obj‏ بیشتر بستر طغیانی رود را شامل می‌شود و این بخش از محدوده در معرض مستقیم 
فرسایش توسط جریان رود می‌باشد. مخصوصاً در محل پیچ مقعر مثاندرهاء به دلیل افزایش سرعت جریان رود و عمل 
زیربری عمل تخریب شدیدتر است و وجود سازندهای آبرفتی حساس از عوامل مهم هستند که رودخانه‌های اصلی 
بیشتر تحت تأثیر بوده است. به علاوه به استناد بازدید میدانی صورت گرفته از منطقه و گفتگو با اهالی بومی 
و مسئولان ذی‌ربط مهم‌ترین شکل solid‏ نامناسب از بستر و حریم رودخانه‌ها و wise‏ رودخانه گیوی‌چای که 
باعث تشدید فرآیند فرسایش و رسوب‌گذاری در این منطقه شده است. شامل؛ توسعهٌ اراضی کشاورزی و EL‏ 
است. از سوی دیگر احداث جاده در منطقه مطالعاتی. با از بین بردن پوشش مرتعی در دامنه‌های پرشیب باعث شده 
در هنگام بارش بهاری و ذوب برف» زمینه جهت فرسایش‌پذیری خاک فراهم گردد. در شکل A‏ نمونه‌هایی از عوامل 


دخیل در افزایش پتانسیل فرسایش‌پذیری خاک در سطح شهرستان گیوی» نمایش داده شده است. 
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شکل (A -٩‏ فرسایش ای مک و jal) ack‏ اد پوشش گیاهی مناسب و محافظ. 8 و ») 

فرسایش کنار رودخانه به علت وجود مواد سست و فرسایش‌پذیر و زیربری توسط رودخانه. (D‏ فرسایش به 
دلیل زیربری و از بین بردن پاشنه دامنه توسط جاده. (E‏ آتش زدن عمدی و خشک کردن پوشش گیاهی منطقه 
جهت تغییر کاربری. ۴) توسعة اراضی کشاورزی و CEL‏ در حریم رودخانه. (G‏ تخریب پوشش ALS‏ طبیعی 


محدوده در اثر احداث سد و همچنین ایجاد راه‌های ارتباطی جدید جهت دسترس به محل سد. 


۶- جمع‌بندی 

تغییرات کاربری اراضیء یکی از عوامل مهم در ایجاد فرسایش خاک است و در سال‌های اخیر, تأثیر متقابل تغییر 
کاربری اراضی و فرسایش خاک به یک نگرانی عمده زیست‌محیطی تبدیل شده است. با توجه به اهمیت موضوع؛ در 
پژوهش حاضر نیز بررسی تغییرات کاربری اراضی و نقش Of‏ در میزان فرسایش خاک در سطح شهرستان گیوی. با 
استفاده از تکنیک‌های دورسنجی و تحلیل چند معیاره MABAC‏ برای دو دوره زمانی ۲۰۰۰ و ۲۰۲۱ موردبررسی 
قرار گرفته است. نتایج حاصل از مطالعهء نشان داد که در سال ۲۰۰۰ کاربری‌های مراتع خوب و متوسط به ترتیب با 
پوشش ۲۶/۶۷ و ۲۲/۱۲ درصد از مساحت شهرستان. بیشترین مقدار مساحت را به خود اخحتصاص دادند و در سال 


سال دوازدهم بررسی اثرات تغییرات کاربری اراضی بر فرسایش خاک .... yY‏ 


۷ درصد از سطح شهرستان می‌باشد. با توجه به نتایج مطالعه در طی باز زمانی مورد مطالعه. مراتع متوسط. 
مراتع خوب پوشش lS‏ متراکم و زراعت آبی کاهش يافته است و در مقابل کاربری‌های مراتع ضعیف. زراعت 
دیم» اراضی پست و هموار (اراضی Cab‏ و مناطق مسکونی افزايش پیدا کرده است. همچنین با توجه به نتایج مطالعه 
می‌توان بیان داشت که فعالیت‌های انسانی» مانند؛ بهره‌برداری غیراصولی از مراتع» تخریب پوششگیاهی در اثر تبدیل 
اراضی مرتعی و جنگلی به زراعی» کشت دیم در اراضی شیب‌دار (ازجمله کشت غلات به‌صورت ديم در اراضی 
شیب‌دار و به‌ویژه کشت در شیب‌های بالای ۱۵ درصد)» شخم خوردن مراتع در شیب‌های تند. عدم رعایت تناوب 
زراعی» استفادة نامناسب از بستر و حریم رودخانه‌ها؛ تغیبر مراتع و BAS‏ آن به مناطق مسکونی و انسان‌ساخت» 
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